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Let us consider the local specification system of Gibbs point pro- 

'^ , cess with inhibition pairwise interaction acting on some Delaunay 

subgraph specifically not containing the edges of Delaunay triangles 

ewith circumscribed circle of radius greater than some fixed positive 
real value R. Even if we think that there exists at least a station- 
rj—k ' ary Gibbs state associated to such system, we do not know yet how 

(-H , to prove it mainly due to some uncontrolled "negative" contribution 

Q ' in the expression of the local energy needed to insert any number of 

^ ' points in some large enough empty region of the space. This is solved 

by introducing some subgraph, called the _R-local Delaunay graph, 
which is a slight but tailored modification of the previous one. This 
►-^ I kind of model does not inherit the local stability property but sat- 

QQ ■ isfies some new extension called i?-local stability. This weakened 

r~^ ' property combined with the local property provides the existence of 

^S| , Gibbs state. 

f— ^ . keywords: Gibbs states, Delaunay triangulation, pairwise interaction, 

\^ ' D.L.R. equations, local specifications, correlation functions. 
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a ■ 1 Introduction 



' O I There exist many different manners to describe Continuum Gibbs models. 

^ ■ One way is using correlation functions [261 127| , another one is rather using 



local specification 25 . One could also investigate integral characteriza- 
tion with Palm distribution |S1|53], empirical measure leading to ergodic 



theorem |24) . variational principle, minimizing the excess free energy den- 
sity linking pressure, entropy and energy density "S". In this framework, 
an important ingredient, particularly useful for the existence of a Gibbs 
measure |25l I26j , is the relative compactness assumption taking different 
form depending on the chosen description (correlation functions, specifica- 
tions and relative entropy...). The local stability is a sufficient assumption 
for the relative compactness assumption. However this property is very 
interesting in several other purposes: convergence of Markov chain Monte 
Carlo (McMC) algorithms to reach equilibrium |14l I17j . stochastic FKG 
domination |11| . uniqueness of a Gibbs state via Kirkood-Salsburg equa- 
tions |22[l26l l5]... One might find interesting to weaken the local stability 
assumption. In the classical framework of pairwise interaction Gibbs point 
process, the natural extension is the well-known superstable assumption 
including hard core, inhibition and Lennard Jones pairwise potential. As 
already discussed for example in W,*¥', this assumption is not well suited for 
nearest neighbours models introduced by Baddeley and M0ller 1 , where 
the neighbourhood relation depends locally on the realization of the pro- 
cess. Kendall et al ^HJ developed models generalizing area interaction 
model called quermass interaction processes and intensively studied simu- 
lation, statistics, Markov properties (see also |16[l^l^) rather in the spirit 
of non stationary Gibbs states. As it is often the case for nearest-neighbour 
continuum models, such a kind of model is not locally stable . 

In jS] we deal with nearest-neighbour continuunr Potts model where 
the repulsion between particles of different type only acts on a Delaunay 



(sub)graph. This work is an adaptation of the Lebowitz and Lieb soft 
core continuum Potts model |19| . In order to exhibit a phase transition 
phenomenon for such a model, we mainly need, on the one hand, to use 
some already known percolation result |15| on the Delaunay graph and, on 
the other hand, to choose some Delaunay subgraph for which the existence 
of the related two-types particles model can be established. The ques- 
tion arises whether one or several of our existing one-type particle nearest 
neighbour model could be adapted to this phase transition problem. The 
spirit of these previous works was to build some Delaunay subgraphs for 
which the local stability holds. Unfortunately, as a direct consequence, 
these resulting subgraphs do not behave locally as the original Delaunay 
graph and do not directly inherit the required percolation property due to 
some change of connectivity when the number of neighbours dramatically 
increases. In order to ensure both the percolation property and the exis- 
tence of the model, the chosen solution was to define the nearest-neighbour 
continuum Potts model as a two- types particles Gibbs point process related 
to some energy function with hard-core component and based on some De- 
launay subgraph specifically not containing the edges of Delaunay triangles 
with circumscribed circle of radius greater than some fixed real value R. 
Let us point out that, without the hard-core assumption, the proof of the 
existence of the one-type particle pairwise Gibbs point process based on 
such a Delaunay subgraph is not yet established even if we additionnaly 
require the nonnegativeness of the interaction function. However, we are 
convinced that this model exists. In other words, this simply means that 



the well-known pairwisc inhibition model on the complete graph is not yet 
adapted to the Delaunay nearest-neighbour framework. The main reason 
is that the local energy needed for the insertion of one point is necessarily 
stable (since obviously nonnegative) for the first one whereas this could 
not happen for the second one due to some "negative" residual edges con- 
tributing in its expression. The present study is an attempt to define a 
first existing version of (pairwise) inhibition nearest-neighbour Gibbs point 
process. We investigate to introduce some new Delaunay subgraph, called 
i?-local Delaunay. In fact, this is a subgraph of the one considered just 
above. It mostly preserves the same expected properties and above all its 
"local" behavior, that is, the same connectivity at small scale as the De- 
launay graph. Its further characteristic is that no "negative" contribution 
of residual edge occurs when inserting any number of points in domain de- 
fined as union of balls of radius R. As a direct consequence, the related 
inhibition nearest-neighbour model inherits some new property, called R- 
local stability, weaker than local stability but enough to upper bound the 
correlation functions of the process by some homogeneous Poisson process 
ones. Combined with the local property intrinsic to such process, this last 
property provides the existence of stationary Gibbs state. 

We may hope that nearest-neighbour continuum models are interesting 
for small temperature (not too small for a classic approach) as an alter- 
native of standard models on regular networks, because it allows more 
degrees of freedom and may be find applications in crystallography. We 
may think of the rigidity and plasticity properties of glasses or the study 



of ferromagnetic fluids or liquid cristals (smectic A,C, neniatic N). See for 
example |3 I12j and references therein. In particular, it seems that the 
place of emptyness is important for the study of equilibrium tension in 
a menbrane [3 I20j . More generally, it is well-known that Voronoi graph 
and regions (rather called Wigner-Seitz grid and Brillouin zone in physics 
framework) take a fondamental place for the understanding of the electrical 
current, waves propagation and phase transitions. 

After giving some notations and preliminaries about the i?-local Delau- 
nay graph in section [21 we introduce in section the inhibition pairwise 
interaction model based on the i?-local Delaunay graph. After introducing 
the definition of the i?- local stability, we establish in section[Slthe existence 
of a Gibbs measure associated to some local specifications family based on 
this related energy function. 

2 The i?-local Delaunay graph 

In this section and for the rest of the paper, R designates some fixed non- 
negative real number. For any given Borel set A C R'^, one denotes by H. 
and r^A the classes of locally finite subsets of points, called configurations 
in this paper, in IR"^ and A respectively. In particular, fi/ denote the sets of 
finite configurations in Q. Moreover, for any set A (not necessarily a Borel 
set) 7^2 (A) designates the set of pairs of points in A. Let B and Bb be the 
set of Borel sets and bounded Borel sets of IR"*. An element (p oiH. could be 
represented as ip ~ y Sxi which is a simple counting Radon measure in 
M"^ (i.e. all the points Xi of M''' are distinct) where VA G S , Sx{A) = 1a (x) 



is the Dirac measure and lyii(.) is the indicator function of a set A. This 
space n is equipped with the vague topology, that is to say the weak topol- 
ogy for Radon measures with respect to the set of continuous functions 
vanishing outside a compact set. J-' is the cr- field spanned by the maps 
if — > '/'(A) , A G Bb, V(p G fl. The corresponding cr-field ^a is similarly 
defined on Q/^ . Furthermore, for any A <E Bb, 

(r!,^) = (r!A,^A) x(r!Ac,^Ac) 

where A^ denotes the complementary of A in _ZR . Let T\ be the reverse 
projection of J^a under the previous identification, so that J-'a is a a-field 

on n. 

A point process ^ on M (respectively <J>a on A) is a random variable 
on Q (respectively on Qa) and is associated to a probability distribution P 
on (rj,^) , (respectively Pa on {Qa,^a))- 

Some configuration ip is said to be in general position when no d + 2 
points lie on the same hypersphere (with no point inside) and no / + 1 
{I — 2, . . . ,d) points lie on some I — 1 dimensional affine subspacc of R"^. 
For any simplex ip (triangle when d = 2) in some configuration 93, one 
denotes by C{ip) the greatest hypersphere circumscribed by ip with no 
point of if inside its interior. The radius and the center (voronoi vertice) of 
such hypersphere are respectively denoted by r{ip) and c(V') • One notices 
that, for any simplex ip, #(C('0) Hep) = d+ 1 holds only if the configuration 
if is in general position. 

Before defining the _R-local Delaunay graph we first need to recall the 
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definition of tlic Dclaunay graph. 

Definition 1 For some ip in Cl in general position, one defines Deld+i{(p) 
by the unique decomposition into simplexes -0 in which the convex hull of 
the hypersphere C('(/') does not contain any point of ip\ip. 
The Delaunay graph is then defined by the set of edges : 

Dehi^p) = U ^2(^). 

iieDela+iiv) 

According to the previous definition, one can assert in the two dimensional 
case that the Delaunay graph is a triangulation whenever the configuration 
(p is in general position. 

Now, we propose to define a subgraph of the Delaunay graph with edges 
that are particularly of length lower than some positive fixed distance R. 
One need first to introduce, for any set A, the set A® R ^ \JxeAB{x, R) 
where B{x, R) = {y G R'^ : ||a; — y|| < R} is the usual ball of radius R and 
||x — y|| designates the Euclidean distance between the points x and y in R'^. 
Moreover, the complement of some set A in IR"^ is denoted by A"^ — R'^\ A. 

Definition 2 For any ip ^ il, one defines : 

1. Tiie i?- vacuum of p : 

^Rif) = U {^(^' ^) • ■^(^' ^) n (^ = 0} = (y^ © i?)' ® R. 

2. Tiie i?-local Delaunay grapii : 

Del^{p>)= Pi DehiipU^) 



where flf{$ii{ip)) is the set of finite configurations in the i?- vacuum 
of ip. 

Thus, one may give further interpretation of this subgraph: any edge of 
the Delaunay graph of ip possibly broken when inserting points in the R- 
vacuum $B.i'i^) of ip does not he in the i?-local Delaunay graph of (p. This 
interpretation leads to another way to define this subgraph. By defining 
the influence region Z^{{x,y}) of any edge {x,y} G De^iv) by: 

ZMx,y}) = n {B{c,r):{x,y}(ZdB{c,r)&ndB{c,r)r\^\{x,y]^{D} 

celR'',r>0 

fl S(c(V'),r(V')) 

i>e-Deid+i(93) : iljD{x.,y] 

one derives another characterization of this subgraph: 

Del^{^) = {{x, y] e DeUip) : Z^{{x, y}) n %R{ip) = 0} . 

Clearly, with respect to this characterization, one may assert some nice 
property: 

Del^'{ip) C Del^^{ip) whenever Ri < R2 

since in this case $ii^{<p) D %i^^{ip) (illustrated in figure^. 

Some other properties of the i?-local Delaunay graph are described in 
the next proposition. 

Proposition 1 For any configurations ip, fi and 1^2, the following prop- 
erties holds : 

1. If %R{if) = then Del^{ip) = Dehiip). 



2. For any R' > R, any point x and any borelien set A , the following 
holds : 

Del^{ip n (A ® RY) C Del^iifi). 



3. When ipi and ip2 are such that d{ipi,ip2) = inf ||a;i — X2II > 



2R, 



Del^iipi U ip2) = Del^iipi) U Del^{ip2). 
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Figure 1: Both figures represent the same configuration (p of some points (or 
vertices). However, the first one is given with the R- local Delaunay graph 
(solid lines) with R — 50 whereas the second one is for R = 100. The other 
edges (dotted lines) are the residual edges of the Delaunay graph. The 
union of the (darker) gray balls corresponds to the quantity i'r{^p) defining 
in some sense the border of the i?-vacuum. The lighter gray part is the rest 
of the i?-vacuum. 

The first property means that, for some configuration <y9, the Delaunay 
graph of if is the same than the i?-local Delaunay graph of ip with R chosen 
great enough. The second one (see figure|2Jl is a key-property since it points 
out that, for this kind of graph, the edges of the graph before the insertion 
of some point inside a region bigger than a ball of radius at least equal 
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Figure 2: The figures represent in solid lines the R- local Delaunay graphs 
{R ~ 200) for the configurations tpO A'^ (on the left) and ip (on the right) 
of some points where A is the Borel set bounded by the square. The set 
of the widest solid edges corresponds to the intersection of both graphs 
Del2{<p n A°) and Del2{f)- In fact, it coincides with the first graph since 
this graph is included in the second one. One may observe that this is 
untrue for the Delaunay graphs represented by solid and dotted edges. 

to R are in the graph after the insertion. The last one asserts that two 
subconfigurations of points separated with a distance greater than 2_R, are 
disconnected in the graph of the whole configuration. 

Another characteristic property of the i?-local Delaunay graph is given 
below. 

Proposition 2 If any {x.y} £ Del2{(p) is such that d{x,y) = \\x — y\\ > 
2R, then {x,y} ^ DelJ^{ip). 

Proof. In this case, for any ip e Deld+i{'^) such that {x,y} C Vj the 
radius r^iji) have to be larger than R ■ 
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More generally, for any {x,y} G DelJ^i^p): we have 

where diam(A) = sup/^ ^ \^/^2 \\zi — Z2II is the diameter of any given Borel 
set A C W^. 

The first property given in the proposition ^ asserts that the i?-local 
Delaunay graph and the Delaunay graph are the same whenever the R- 
vacuum is an empty set. Some kind of generalization of this property is 
given below by describing the local behavior of the i?-local Delaunay graph 
observed on some particular region of R'^ related to the value of R. 

Proposition 3 For any configuration ip and any Borel set A such that : 

which is satisfied whenever (A ® i?) n 9r{(p) — 0, one may assert that : 

Clearly, the i?-local Delaunay graph coincides with the Delaunay graph in 
region of R'' with a concentration of points large enough with respect to 
the value of R. 

Now, for pratical purposes, we attempt to give an equivalent definition 
of the i?- local Delaunay graph by introducing a subregion of the i?- vacuum. 
Indeed, in the current form the i?-local Delaunay graph is not easily com- 
putable. We then introduce a special set of edges : 

£Riip)= U \J{^\{x}:3ze^\{x},rii^\{x})U{z})<R}. 

ijeDeld+i{<p):r{Tp)>R xGip 
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For any ^ G £ii{(p), let us denote by cr{^) the unique point such that : 
(Vx e e, \\x - ciiiOW = R) and B{cniO,R) n{ip\0 = 9. 

We propose to define some sort of border of the _R- vacuum 0_r(<p) of some 
configuration ip by introducing the subset Vr{(p) of 9r{(p) : 

Vniv)^ U B{cr{0,R). 

Consequently, the i?-vacuum can be decomposed into two parts (not nec- 
essarily disjoint) : 

^R{ip)^VR{ip)ul U i?(c(V),r(V')) 

\-4>eDela+i{v)-rW>R 

This is illustrated by the previous figures in the two-dimcnsionnal case. 

By applying the following proposition, it is almost easy to propose an 
algorithm in order to compute the i?-local Delaunay graph for some con- 
figuration (p. 

Proposition 4 The residuals edges in the difference between the Delaunay 
graph and the i?-local Delaunay graph can be decomposed into the union 
of two disjoint sets : 

Dehi^) \ Del^i^p) = nes^{Lp) U 7^esf (93) 

where 

Teesf (^) = U {C e 7^2 (V') : r(V) > R} 

and 

T^es^iifi) == {e e Dehiip) \ Ues^iip) : Z^ (C) n Vr (ifi) + 0} 
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Remark 1 The same idea can be applied to some classical subgraphs of 
the Delaunay graph like the Gabriel and the Relative Neighbours graphs. 
Recall that these subgraphs are respectively defined as follows: 

G(^) = \{x, y} e Deh{^) : B{^, t^Jil) n^\{x,y} = 



and 

RNG{^) = {{x, y} e Deh{ip) : B{x, \\x - y\\) n B{y, \\x - y\\) n ^ \ {x, y} = 0} . 

Indeed, we have just to define the infiuence regions Z^({x, y}) and Z^^^{{x, y}) 
of any edge {x, y} of each graph respectively by : 

and 

Z^''''{{x,y}) ^ B{x, \\x -y\\)n B{y, \\x - y\\). 

One may derive one another characterization of these new i?-local sub- 
graphs by asserting : 

G«(^) = {{x, y} G GM : Z^{{x,y}) n 0k(^) - 0} . 

and 

RNG^i^) = {{x, y} e RNG{^) : Z^^«({x, y}) n 0;,(^) - 0} . 

3 Inhibition interaction model on the i?-local 
Delaunay graph 

We first introduce the definition of the energy function induced by the 
i?-local delaunay graph. 
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Definition 3 Given any fixed i? > 0, one defines the i?-energy of some 
finite configuration Lp by : 

where (\> is some upper bounded nonnegative pairwise interaction function. 

The previous model could be easily extended by adding interaction terms 
function of all order. In the rest of this paper, we only deal with the 
pairwise interaction but all the results remains valid for these extensions 
whenever the interaction functions of all order are nonnegative and upper 
bounded. 

Remark 2 Another similar model is the one with pairwise interaction be- 
tween Voronoi vertices. In fact, each influence region of Delaunay edge 
is the intersection between two Delaunay disks for edge inside the convex 
hull of points and just one Delaunay disk for edges belonging to the convex 
hull. The local Delaunay graph then take into account interaction between 
Voronoi vertices and some points characterizing the R-vacuum region. By 
denoting V((^) = {c('0)}v,GDei3((ij) '^^ '''^^ of the Voronoi vertices we define 

Vi'M = {{c(V'i),c(V'2)}e(vM)':#(V'in7^2) = 2, r(^i) < i?, r(V'2)<i?, 

i?(c(Vi), r(V'i)) n BWV-z), r(^2)) n B{cR{i), R) = 0, V^ £ £r{^>)\. 
The following finite energy is of the same kind of the i?-energy: 

yR{^)^ E 0({c(V'i),c(V2)}) 

{c(^l),c(^2)}eV2«(¥') 
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where (f> is some upper bounded nonnegative interaction function between 
Voronoi vertices. 

As usually, the mutual energy and the conditional energy between two 
configurations ip and tp are respectively defined by : 

W^ {(p, tp) ^V^ifUtp)- V^ [if) - V^ (-0) 

and 

Due to the properties of the i?-local Delaunay graph, we have the fol- 
lowing result. 

Proposition 5 f . For any set A, one has : 

V'^ ('/'Aefl|'/'(Aefl)0 == y\ 0({a;,y}) (2) 



{a;,y}eDe/«(vJ)\De;«(<^(AeH)<=) 



and 



{a;,j/}e_De/«(<^)\(De;«(ip(Affifl)<=)UDe/«(ipA®B) 

2. ff Ai and A2 are two Borel sets such that d(Ai,A2) > 2_R, the 
following holds : 

3. Finite range property: for any bounded Borel set A, 

V" (</^a|</'A=) - V" (<^A|<^Ae6fl\A) • (3) 
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Proof. In fact because of the translation invariance property of the local 
energy, it is sufficient to prove that for any <y9 G f2, 

y«(o|(^) = y^(o|^6i?)- 

where, for shortness, one denotes (pgi?, = fBio.eR) ■ Roughly speaking, let 
us show that 

Del^i^)ADel^{^ U {0}) = Del^i^eR)^Del^i^eR U {0}) 

where A denotes the symmetric difference operator (i.e. for two sets A 
and B, AAB = {A\B)Ll{B\ A)). 

Clearly as one remarks after proposition 4, the diameter of the influence 
region for any edge of the _R-local Delaunay graph is bounded by 2R. 

So one has for any ^ — {x,y} G Del2(Lp)A.Delf^{(f U {0}) (resp. ^ = 
{x, y} G Del§i^eR)^Del§{{ipeR U {0}))): 

where v can be chosen as (p or ipU {0} (resp. (pen or ip^ji U {0}). Moreover 
if one notices that 

0k((^) n B(0, AR) =. 9r{^6r) n B{0, AR) 

and 

0fl(^ U {0}) n 5(0, 4R) - 9R{ipeR U {0})) n 5(0, 4R) 

it follows by the definition of the _R-local Delaunay graph that 
Del§{ip)ADel^{ipU{0}) = Del^{ip)ADel^{ipU {0}) nP2{B{0,4:R)) 
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= Del^{^^B)^Del^{^eB, U {0}) n ^2(^(0, 4i?)) 
= Del^{^^R)J^Del^{^eR U {0}) 

and the proof is complete H 

Defined as a subgraph of the Delaunay graph, the i?-local Delaunay 
grapli is linear in the planar case and then the i?-energy V^ inherits of the 
stability property. In the higher dimensional case, stability occurs since the 
interaction function is assumed to be nonnegative. 

Furthermore, one can notice that the i?-energy is not locally stable due 
to its local behavior as the Delaunay graph on each region of space with 
high enough density of points. This present work is then really different 
from the previous ones 01 |3j where the goal was to build some subgraphs 
of the Delaunay graph providing local stability. 

4 The R-\oca\. stability 

However, one may assert some new property for the i?-energy V^ which is 
an extension of the local stability. 

Definition 4 For some nonnegative real number i?, an energy function V 
is said to be i?-locally stable if there exists K >Q such that for any finite 
configurations 93, ipi and ip2, and any subset A satisfying A 3 •^,^i C A(BR 
and ¥^2 C (A © Ry-. 

V{ipLI ipi\ip2) > -K#{ip U ipi) (4) 
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One may arrange the i?-local stability as some property between the stabil- 
ity (acting globally) and the local stability. Indeed, on the one hand when 
R vanishes and A = (p, the i?-local stability is similar to the local stability 
and on the other hand when A = R'' (i.e. (p2 — 0) the i?-local stability 
coincides with the (global) stability. 

Proposition 6 The i?-energy V^ with nonnegative interaction function 
is i?-locally stable. 

Proof. This is a direct consequence of the property Q with K = M 

This property will play later some key-role in the proof of the non 
emptyness of stationary Gibbs state based on the i?-energy. 

Given any outside configuration (p°, Ruelle [201 has introduced the fol- 
lowing quantity: 



Zi (^°) /a\a 
where for any measurable function / : 



® dipf{ip)^'^— dxi .. .dXnfixi,.. .,Xn)- 

Ja n. Ja" " » ' 

As a particular case, A = 0, one can derive the correlation function p^ {ip\ip°) 
PlAiflv"), satisfying 

An interesting well-known property is then to prove that this correlation 
function are upper bounded by the correlation function of some Poisson 
process, that is, of the form ^'^'^. 

18 



Proposition 7 If some energy function V is i?-locally stable tlien 

Pa {v\^°) < C*^ wliere ^ = ze^e^^^'I^^O'^)!. 
Proof. For any A and A D A © i?, 



■^A ('^°) /An(A©_R) /An(Aefl) 

^ \z( f diPie'^*'!'' i dV2e~^('^^l'^°) (by i?-local stability) 

^A 1'/'°) JAn(Ae_R) JAn(A©fl)<= 

< f-c-^^^'^ ze-^|An(Aefl)| ^An(Aefl)'= (y^°) 

- ^ ^ zx(<^°) 

< (ze^)*^ ,--|An(Ae«)| (^j,,, Zl,(^^^). (^°) < Z^ (^°) ) 
Finally, applying this result when A = </?, this leads to 

since |A n (A © i?)| is then upper bounded by #(/? x |i3(0, R)\ ■ 

5 Existence of a Gibbs measure based on the 
i?-energy 

At this stage, everything that one needs in order to prove the existence of a 
Gibbs measure related to the energy function y^, was already introduced. 
We then first recall the definition of local specifications based on the R- 
energy. 

Definition 5 Given any outside configuration (p°, the following family of 
measures 
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n = {n\{- , uj)} \^i3^ on {^,T) is a system of local specifications : 

^A(^) J A 

where the partition function is given by Z^ {ip°) = (p dip e~ 



VMvla) 



In order to prove existence of Gibbs state related to this system of local 
specifications, the following probabilties HaI-PaI'^") for any bounded Borel 
sets A C A, and any Fa G Ta, have to be controlled uniformly on A D A 
and (p°. By denoting F^"^ the projection of Fa onto Ta, 



n»(F.b") . £.^.,..m(^£, 



^ d^lpi..{^) -^ (6) 

where Pa a if\v°)j defined in Q, play an important role in the expression 
of the Radon- Nikodym of the local specification nA(-|<p°) with respect to 
some Poisson process in A. 

Consequently, by combining the result of the proposition [3 with K = 
and the equation © one derives for any bounded Borel sets A C A, and 
any Fa G J- a that, 

Pa,a (V'l'/'") 



TIa{Fa\^°) = (P diplp,^..{ip)- ^^^ 



< j,^ d^lp,^.{^)- ^^^ 



3(0, H) 

d^lj.,oc(^) (7) 



/A 

In particular, for the event Fa ~ [<I>a(A) > to] one obtains 



+ 00 fc « +00 , z|B(0-R)||A 

pa([$a(a) > to]|^°) = ^ i^ y^^ rfa:i . . . dx.pi^A ({^1, • ■ • , ^™}i^°) < E — — ^T"^ 



k—m k=m 
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In some sense, this means that Pa is "dominated" by the non-normaUzed 
Poisson process with intensity 2;e^l'^'-'^'''^'l|A|. 

In 0] , we proposed some simpler sufficient conditions based on the local 
energy in order to satisfy the Preston's theorem (|5S] theorem 4.3 p. 58) 
assumptions very useful for proving the existence of a stationary Gibbs 
state. 

(LS) Local Stability: there exists some constant K > : 

V{0\ip)>~K, Vipen (8) 

(Q) Quasilocality: for any bounded Borel set such that G A : 

\V{0\p) - F(0|<^a)| < £(d(0, A^)),V<^ G n (9) 

where e is a nonnegative decreasing function which vanishes asymp- 
totically and d{x, B) = mind(a;, y) is the Euclidean distance between 
a point X and a Borel set B. 

The first assumption is not satisfied by our model based on the R-energy. 
Fortunately, one may replace it by the i?-local stability assumption: 

(i?-LS) i?-local stability: there exists some real value R > Q such that 
V is -R-locally stable. 

but also by the more general one based on the correlation function: 

(UC) Upperbound of correlation function: there exists some real value 
5 such that p%{Lp\if'') < £.*'^. 
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Proposition 8 By assuming that 11 is a system of translation invariant 
local specifications based on some energy function satisfying (UC) (implied 
by (R-JjS)) and (Q), the set ^o(n) of stationary Gibbs measures is non 
empty. 

Proof. The proof of this result is very similar to the one proposed in ^. 
The only difference is that condition (3.7) of theorem 4.3 (in [21] p. 58) is 
satisfied as a direct consequence of equation {T)) ■ 

Consequently, we may assert the main result of this paper. 

Theorem 9 The set ^o(n''^) of stationary Gibbs measures associated to 
the system II-'^ of translation invariant local specifications based on the 
i?-energy V^ defined in (Q is non empty. 

Proof. V^ is i?-locally stable (see proposition O and satisfies the finite 
range property (of proposition [SI H 

We finally end this section by some concluding remarks. 

Remark 3 In the plane, the maximum number of Voronoi vertices (or 
equivalently, the number of Delaunay triangles) is an upper bound for the 
number of holes (or more precisely, the Eulcr characteristic) generated by 
the Quermass-interaction model studied in |18j . Thus, there is a strong 
link between our model and the quermass-interaction model in the planar 
case when the grains are disks of fixed radius. The Quermass-interaction 
model ,18 and nearest neighbours models [l defined using the Delaunay 
graph raised problems of stability in dimension greater than two whereas 
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models presented here works in any dimension. 

Remark 4 The assumption of relative compactness we discussed here for 
some particular nearest neighbours models is useful when we use corre- 
lation functions or local specifications. This assumption appears also in 
the modern large deviation theory and for having the sub levels of the 
specific entropy sequentially compact and existence of an accumulation 

point nnnnmni. 

Remark 5 In ^, in order to study phase transition, we introduce the 
nearest-neighbour continuum Potts model where the soft repulsion between 
particles acts on a graph defined by the Delaunay edges. In order to prove 
the existence of this model, we simply add (in the spirit of [3]) an hard-core 
component acting on all particles independently of their type. However, 
this result is still true without this assumption using the -R-local Delaunay. 

Remark 6 The stability of the finite energy and the temperedness of the 
mutual energy (see |26j p. 32), implied by finite range property, provide 
results of |2S| (p. 41-58) concerning the existence of the pressure with free 
boundary condition and thermodynamic limit of microcanonical, canonical 
and grand canonical ensembles. 
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